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Abstract 


The  broadband  scat  tering  characteristics  of  a  target  may  be  used  to  distinguish  echos 
from  a  target  of  interest  from  those  of  clutter.  This  type  of  classification  is  of  interest 
in  many  sonar  applications:  ASW,  torpedo,  surface  mines,  and  diver  detection.  In 
a  shallow  water  situation  the  original  sonar  pulse  and  the  echo  from  the  target  will 
consist  of  a  sequence  of  pulses  corresponding  to  the  various  combinations  of  incident 
and  backseat tered  multipath  arrivals.  For  a  very  shallow  target  or  one  close  to  the 
seabed,  it  may  be  difficult  to  isolate  some  of  these  multipath  arrivals.  It  is  expected 
that  these  waveguide  effects  could  have  a  significant  effect  upon  some  classifiers.  In 
this  paper,  we  discuss  the  modelling  of  scattering  of  sound  from  an  elastic  sphere 
(shelled,  solid,  water-filled,  etc)  which  is  in  a  waveguide  with  an  upper  surface  and  a 
lower  seabed.  An  exact  propagation  solution  is  implemented  and  used  to  benchmark 
a  ray-expansion  solution.  The  scattered  signals  for  a  variety  of  different  sonar  pulses 
and  different  positions  of  a  spherical  target  in  the  waveguide  are  shown. 

Resume 


On  peut  se  baser  sur  les  caracteristiques  de  diffusion  a  large  bande  d’une  cible  pour 
faire  la  distinction  entre  les  echos  d’une  cible  d’interet  et  le  clutter.  Ce  type  de  clas¬ 
sification  est  utile  dans  bon  nombre  duplications  sonar  :  guerre  anti-sous-inarine, 
defense  contre  les  torpilles,  chasse  aux  mines  de  surface  et  detection  des  plongeurs.  En 
eau  peu  profonde,  Fimpulsion  sonar  initiate  et  l’echo  produit  par  la  cible  consistent 
en  une  sequence  d’impulsions  correspondant  aux  diverses  combinaisons  de  signaux 
incidents  et  de  signaux  retrodiffuses  regus  par  trajets  multiples.  Dans  le  cas  d’une 
cible  a  tres  faible  profondeur  ou  d’une  cible  a  proximite  du  fond  marin,  il  peut  etre 
difficile  d’isoler  certains  de  ces  signaux  regus  par  trajets  multiples.  On  prevoit  que 
les  effets  de  guide  d’ondes  pourraient  avoir  une  incidence  appreciable  sur  certains 
systemes  de  classification.  Le  present  document  traite  de  la  modelisation  de  la  diffu¬ 
sion  du  son  par  une  sphere  elastique  (a  coquille,  pleine,  remplie  d’eau,  etc.)  comprise 
dans  un  guide  d’ondes  avec  une  surface  superieure  et  une  surface  inferieure  constitute 
du  fond  marin.  Une  solution  de  propagation  exacte  est  mise  en  oeuvre  et  utilisee  pour 
servir  de  reference  a  une  solution  de  prolongement  des  rayons.  Les  signaux  diffuses 
pour  differentes  impulsions  sonar  et  differentes  positions  d’une  cible  spherique  dans 
le  guide  d’ondes  sont  presentes. 
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Executive  summary 


Modelling  broadband  scattering  from  shelled  spheres  in  a 
waveguide 

John  A.  Fawcett;  DRDC  Atlantic  TM  2007-270;  Defence  R&D  Canada  -  Atlantic; 
October  2007. 

Background 

The  detection  and  classification  of  a  sonar  echo  is  a  problem  of  much  interest  for  many 
sonar  applications:  ASW,  torpedo  defence,  minehunting  and  diver  detection.  The 
ability  to  model  a  broadband  signal  scattered  from  a  target  in  the  oceanic  waveguide 
is  important  in  simulating  data  for  sonar  detection  and  classification  studies.  In  this 
report,  we  consider  a  spherical  scattering  object  located  within  a  constant  sound 
speed  waveguide  with  an  upper  surface  and  lower  seabed. 

Principal  results 

It  was  found  that  it  was  possible  to  accurately  and  efficiently  model  the  scattering 
from  a  spherical,  shelled  structure  in  a  waveguide  for  a  bandwidth  of  80  Khz  (or  more). 
The  effects  of  the  waveguide  on  the  resultant  echo  structure  can  be  significant. 

Significance  of  results 

It  is  shown  that  in  interpreting  the  sonar  echos  from  objects  in  a  waveguide,  that 
the  boundary  interactions  of  the  incident  and  scattered  fields  must  be  considered. 
The  model  of  this  report  provides  a  way  of  accurately  computing  the  broadband 
response  of  elastic,  shelled  spheres  in  a  shallow  water  waveguide.  The  spectral  content 
of  the  echo  from  a  target  can  vary  significantly  with  its  depth  and  range  in  the 
waveguide  or,  in  the  time  domain,  the  coherent  combination  of  the  echos  from  the 
various  multipaths,  can  change  the  received  pulse  significantly  from  that  which  would 
be  received  in  free  space.  This  could  lead  to  significant  effects  in  terms  of  target 
detectability  and  could  affect  the  classification  of  the  object  in  the  case  that  the 
classification  was  based  upon  the  free  space  response  of  the  target. 


Future  work 

In  future  work,  we  would  like  to  use  the  model  of  this  report  to  generate  data  for 
detection  and  classification  studies.  Some  simple  experiments  could  establish  the 
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validity  of  the  model.  The  model  could  be  extended  to  more  general  sound  speed 
profiles  by  interfacing  it  with  a  ray-tracing  code. 
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Modelling  broadband  scattering  from  shelled  spheres  in  a 
waveguide 

John  A.  Fawcett;  DRDC  Atlantic  TM  2007-270;  R  &  D  pour  la  defense  Canada  - 
Atlantique;  octobre  2007. 

Introduction 

La  detection  et  la  classification  d'un  echo  sonar  constituent  un  probleme  qui  preseute 
beaucoup  d’interet  dans  bon  nombre  duplications  sonar  :  guerre  anti-sous-marine, 
defense  contre  les  torpilles,  chasse  aux  mines  et  detection  des  plongeurs.  La  capacite 
de  modeliser  un  signal  a  large  bande  diffuse  par  une  cible  dans  le  guide  d’ondes 
oceanique  est  importante  pour  la  simulation  des  donnees  destinees  aux  etudes  de 
detection  sonar  et  de  classification.  Dans  le  present  rapport,  nous  considerons  un 
objet  diffusant  splierique  compris  dans  un  guide  d’ondes  a  vitesse  du  son  constante 
avec  une  surface  superieure  et  une  surface  inferieure  constitute  du  fond  marin. 


Resultats 

On  a  constate  qu’il  etait  possible  de  modeliser  de  nianiere  precise  et  efficace  la  dif¬ 
fusion  par  une  structure  spherique  a  coquille  comprise  dans  un  guide  d’ondes  pour 
une  largeur  de  bande  de  80  kHz  (ou  plus).  Le  guide  d’ondes  peut  avoir  une  incidence 
appreciable  sur  la  structure  de  l’echo  resultant. 


Portee 

On  montre  que  dans  Interpretation  des  echos  sonar  produits  par  des  objets  dans  un 
guide  d’ondes,  les  interactions  aux  limites  des  champs  incidents  et  diffuses  doivent 
etre  prises  en  compte.  Le  modele  du  present  rapport  permet  de  determiner  de  maniere 
precise  la  reponse  a  large  bande  de  spheres  a  coquille  elastiques  dans  un  guide  d’ondes 
en  eau  peu  profonde.  Le  contenu  spectral  de  l’echo  produit  par  une  cible  peut  varier  de 
fagon  appreciable  en  fonction  de  sa  profondeur  et  de  sa  distance  dans  le  guide  d’ondes 
ou,  dans  le  domaine  temporel,  la  combinaison  coherente  des  echos  regus  par  divers 
trajets  multiples  peut  donner  une  impulsion  grandemeut  differente  de  ^impulsion  qui 
aurait  ete  regue  en  espace  fibre.  Cette  situation  peut  avoir  des  effets  marques  sur  la 
detectabilite  des  cibles  et  influer  sur  la  classification  de  l’objet  lorsque  la  classification 
est  basee  sur  la  reponse  en  espace  fibre  de  la  cible. 
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Recherches  futures 


Dans  les  recherches  futures,  nous  aimerions  utiliser  le  modele  du  present  rapport  pour 
generer  des  donnees  destinees  aux  etudes  de  detection  et  de  classification.  Certaines 
experiences  simples  pourraient  valider  le  modele.  Le  modele  pourrait  etre  etendu  a 
des  profils  de  vitesse  du  son  plus  generaux  par  son  interfagage  avec  un  code  de  trace 
de  rayon. 
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1  INTRODUCTION 


Over  the  last  few  years,  there  has  been  much  interest  in  the  discrimination  of  objects 
of  interest  from  clutter  on  the  basis  of  the  spectral  features  of  the  received  echo. 
This  has  been  investigated  for  buried  and  proud  mine  and  general  target /clutter 
classification  [1-4].  Unfortunately,  a  target’s  location  within  the  ocean  waveguide 
will  also  affect  its  echo  characteristics.  For  example,  in  Ref.  5  it  was  shown  how  the 
echo  from  an  aluminum  sphere  on  a  seabed,  was  significantly  different  in  amplitude 
and  phase,  from  that  expected  in  free-space.  Thus,  an  understanding  of  the  target 
scattering  process  in  a  waveguide  is  important  in  understanding  target  detection  and 
classification  in  a  shallow  water  scenario. 

In  this  report,  the  scattering  solution  for  a  sphere  in  a  free  space  (i.e.,  just  surrounded 
by  a  single  homogeneous  fluid)  is  described.  The  sphere  can  have  an  elastic  shell  and 
a  different  interior  elastic  fill.  Thus,  the  model  of  the  sphere  is  quite  general  and 
can  model  a  possible  target  (for  example,  a  steel  shell  surrounding  an  evacuated  or 
filled  interior)  and  clutter,  a  solid  object,  representing,  perhaps,  a  rock.  Some  sample 
spectral  and  pulse  computations  are  presented.  This  model  is  very  efficient  and  can  be 
used  to  model  the  spectral  response  of  spherical  targets  up  to,  for  example,  hundreds 
of  kilohertz  [6]. 

The  exact  solution  and  approximations  for  target  scattering  in  a  waveguide  have  been 
discussed  by  a  variety  of  authors  [7-10].  Similar  expression  to  the  one  in  this  report 
for  the  wavenumber  approach  can  be  found  in  [7-9].  However,  this  method  begins  to 
become  computationally  expensive  for  high  frequencies  or  long  ranges.  In  this  report, 
we  use  a  multipath  ray  expansion  to  describe  the  incident  and  scattered  fields.  A  ray 
type  approach  was  also  used  in  Ref.  11  to  model  scattering  from  a  rigid  sphere  in 
an  oceanic  waveguide  and  is  used,  in  general,  in  many  sonar  performance  models.  In 
the  present  work,  a  systematic  multipath  expansion  is  used  for  the  computation  of 
the  incident  and  scattered  fields.  The  targets  are  taken  to  be  general  elastic-shelled 
spheres,  possibly  with  an  interior  fill.  The  same  general  technique  could  be  used  with 
other  target  shapes,  but  the  use  of  spheres  allows  for  the  scattering  problem  to  be 
solved  analytically  and  the  symmetries  of  the  spherical  geometry  make  some  of  the 
computations  easier. 

In  the  numerical  examples,  some  example  spectral  and  pulse  computations  are  consid¬ 
ered  for  a  steel-shelled  sphere  in  free  space.  Then  the  sphere  is  considered  in  a  shallow 
waveguide  with  a  bottom,  penetrable  seabed  and  a  top  pressure  release  surface.  A 
point  source  is  considered  at  a  specified  depth.  The  resulting  scattered  pressure  fields 
as  computed  by  both  the  wavenumber  integral  and  the  multipath  ray  approaches  are 
shown  as  a  function  of  the  sphere  depth.  In  the  case  of  the  wavenumber  integral  ap¬ 
proach,  the  effect  of  including  rescattering  terms  is  also  investigated.  The  ray  method 
is  then  used  to  compute  the  received  echo  time  series  for  various  ranges  and  depths 
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2  Theory 

2.1  Scattering  from  a  sphere  in  free  space 

The  pressure  field  incident  upon  a  sphere  in  free  space,  due  to  a  point  source  at 
( Rs ,  </><?,  Os),  can  be  expressed  in  terms  of  spherical  Hankel/Bessel  and  Legendre  Func¬ 
tions  as 


OO  00 

Pulc{R,(f),9 )  =  ik  mhn(kRs)jn{kR) 

m— 0  n—rn 

x  P^(cos<j>s)Pn  {cos ^>) cos m(9  -  6 s)  0) 

where  k  =  2nf/c  and  em  =  2  for  m  >  0  and  eTn  =1  for  m  =  0.  In  this  report  we  use 
the  normalization  that 


P?(x)  =  v/2nTI 


‘{n  -  m)!  pm(  \ 
(n  +  m)!  "  [  ’ 


(2) 


where  P”L(x)  are  the  standard  Legendre  Functions.  The  expansion  of  Eq.(l)is  based 
upon  the  origin  of  the  coordinate  system  being  at  the  centre  of  the  sphere.  The 
variable  R  is  the  three-dimensional  range.  The  variable  (f)  is  the  polar  angle  measured 
off  the  vertical  angle  and  9  is  the  azimuthal  angle  around  the  vertical  axis.  The 
subscript  S  denotes  the  coordinates  of  the  source  as  measured  from  the  origin.  In  the 
case  that  4 }s  is  taken  to  be  zero  (i.e.,  the  source  is  taken  to  lie  upon  the  axis  passing 
through  the  sphere)  then  this  expression  simplifies  to 

OO 

pinc(R,<P)  =  ik  \/(2n  +  i)hn{kRs)jn{kR)Pn{cos <t>).  (3) 

n— 0 

The  field  scattered  by  a  sphere  is  then  given  by,  in  the  case  of  Eq.(l) 

OO  OO 

PSC(R,  <M)  =  E  E  «"PnmC‘(cOS  4>)hn{kR)  (4) 

rn~ 0  n— m 


or,  for  Eq.(3). 


pSC(R,(t>)  =  J2anPnChn(kR)Pn(cOS<t>). 

n~  0 


(5) 


In  these  equations  the  terms  p"lr^  and  p\[lc  denote  the  incident  coefficients  of 
jn(kR)p™(cos  4>)  as  defined  by  the  expressions  of  Eqs.(l)  and  (3).  The  scattering 
process  of  the  sphere  is  characterized  by  the  scattering  coefficients  an.  For  example, 
in  the  case  of  a  rigid  sphere, 


h'n{kRo)  ' 


(6) 
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In  the  case  of  an  elastic  shelled  sphere  with  possibly  an  internal  elastic  fill,  a  vector  of 
coefficients  are  obtained  by  solving  a  7  x  7  system  of  equations  for  each  order  n .  There 
are  the  unknown  coefficients  of  the  spherical  Hankel  functions  in  the  surrounding  fluid. 
The  compressional  and  shear  potentials  in  the  shell  each  have  coefficients  for  spherical 
Hankel  and  Bessel  functions.  In  the  interior,  there  are  unknown  coefficients  for  the 
spherical  Bessel  functions  for  the  compressional  and  shear  potentials.  There  are  no 
spherical  Hankel  functions  here  as  they  are  singular  at  the  origin.  The  continuity  of 
displacements  and  stresses  at  the  interfaces  yield  the  necessary  equations  to  determine 
all  the  coefficients  at  each  order  n.  In  particular,  we  are  interested  in  the  coefficients 
of  the  spherical  Hankel  functions  in  the  fluid  surrounding  the  sphere.  Once  these 
coefficients  have  been  obtained  then  the  scattered  field  can  be  computed  anywhere 
in  the  surrounding  free  space. 

However,  when  the  sphere  is  in  a  waveguide,  the  situation  is  more  complicated.  The 
field  incident  upon  the  sphere  includes  the  propagation  effects  of  the  energy  in  the 
waveguide.  The  field  scattered  by  the  sphere  also  interacts  with  the  waveguide  before 
arriving  at  the  receiver.  In  addition,  some  of  the  scattered  energy  will  reflect  off  the 
boundaries  and  be  reincident  upon  the  sphere  and  this  will  lead  to  a  sequence  of 
interactions  between  the  sphere  and  the  boundaries. 

2.2  Scattering  from  a  sphere  in  a  waveguide  -  an  exact 
wavenumber  integral  expression 

In  this  section,  we  derive  an  expression  for  the  scattering  from  a  sphere  in  a  waveguide. 
This  problem  has  been  addressed  in  various  forms  by  previous  authors  [5,7-10].  The 
derivation  below  leads  to  the  various  wavenumber  integrals  which  we  implemented  nu¬ 
merically.  In  order  to  derive  the  scattering  from  a  sphere  in  a  waveguide,  we  consider 
as  a  starting  point,  Eq.(l),  for  the  pressure  field  scattered  by  the  sphere.  We  take 
the  polar  angle,  0,  to  be  measured  off  the  vertical  axis  z.  The  term  hn(k R)P™ (cos  <f>) 
can  be  written  for  <2  <  0  as 

hn(kR)P?( cos  0)  =  r-m  ^ <2j^—P™(l)dq  (7) 

and  for  z  >  0. 

hn(kR)P?( cos^)  =  in~m  (8) 

where  k  =  2i xf/c  and  7  =  \Jk2  —  q2.  These  two  expressions  represent  downgoing 
and  upgoing  wavefields  respectively.  Because  of  the  interactions  with  the  top  and 
bottom  boundaries  of  the  waveguide,  these  fields,  for  each  value  of  the  wavenumber 
parameter  q  will  be  modified  by  a  series  of  reflections  off  the  top  and  bottom  of  the 
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waveguide.  There  will  be  4  wavenumber  integrals  expressing  the  propagation  of  the 
(nm)  Spherical  Harmonic, 


Ti = «^_m_Pir(7A.)exp(!7(2D + 2rre))(,,,  (9) 

T2  _  p(.7(2>,  -  *.>»,  do) 

T3  =  jT  q-^  T  —  ^)(!,27„  C(-7/t)  ex[,(i7(2ft  -  2D  -  2,„) )</,,  (11) 

d 

7-1  =  r~"  r  ^7  0  1  ^ "PPPl*  +  (12) 


Here  /i  is  the  depth  of  the  waveguide,  D  is  the  distance  of  the  sphere  from  the 
seabed,  and  zrcc  is  the  vertical  coordinate  of  the  receiver,  relative  to  the  sphere 
centre.  The  term  /2(g)  is  the  reflection  coefficient  for  the  bottom  seabed  and  for  a 
simple  homogeneous  seabed  has  the  form, 


m  =  1ml zfnA. 

(mi  +  Pi  72) 


(13) 


with  7j  =  y'u^/rj  —  (/2  and  7j  =  y/u2 /c\  —  q 2.  In  the  propagation  of  the  scattered 
field,  each  of  the  above  terms  is  multiplied  by  arlptnc(nm),  the  incident  field  coefficient 
for  the  (nm)  Spherical  Harmonic  multiplied  by  the  sphere  scattering  coefficient  an. 
For  the  rn th  azimuthal  order,  these  propagation  terms  are  summed  with  respect  to 
the  index  n.  Thus,  for  example,  the  summation  of  T\  results  in  the  propagation 
integral 


T \ 

where 


l 


00  qim+lFtn(q)Jm(qr)  R(q) 


k  1  +  R(q)et2lh 

v  Pm  (  T'Y  /  fr\ 

Fm (q)  =  ^  a„— - - —  cos(mn)ptnc(nm) 


Knh/k)  exp(*7(2D  +  zTrc))dq ,  (14) 


7 


(15) 


In  Eq.(15)  the  negative  sign  is  used  with  the  Legendre  Function,  if  the  sign  is  positive 
in  Eqs.(9)-(12)  and  vise  versa.  We  chose  to  write  the  expression  for  Fm  in  this  manner 
as  it  is  the  Scattering  Function  for  the  sphere,  for  the  in  th  azimuthal  order,  in  the 
polar  direction  (f>  =  7 /k,Q  =  7 r. 


DRDC  Atlantic  TM  2007-270 


5 


The  first  of  the  4  expressions,  Eq.(9)  represents  a  reflection  off-  the  seabed  of  the 
scattered  field  and  then  1,2,...  reflections  off  the  water/air  and  seabed  interfaces  before 
arriving  at  the  receiver.  The  second  integral,  Eq.(10),  is  the  same  except  that  there  is 
an  addit  ional  reflection  off  the  top  pressure  release  surface  before  the  scattered  energy 
arrives  at  the  receiver.  The  third  and  fourth  terms  represent  energy  initially  scattered 
upwards  from  the  sphere  with  a  sequence  of  reflections  off  the  top  and  bottom  surfaces 
before  arriving  at  the  receiver  with  an  additional  final  reflection  off  the  bottom  seabed 
in  the  fourth  expression.  There  is  also  the  direct  term,  representing  the  free  space 
scattering  term  which  does  not  undergo  any  interaction  with  the  boundaries  before 
reaching  the  receiver  and  for  this  we  could  use  a  wavenumber  integral.  However,  we 
simply  use  an  expression  of  the  form  of  Eq.(5)  for  this  term.  The  final  expression  for 
the  backscattered  field  is  then  given  by 

psc(r ,  zrec)  =  T\  +  Tjj  +  T3  +  T4  +  To  (16) 

where  To  is  the  direct  propagation  of  the  scattered  field  to  the  receiver.  This  direct 
term  can  be  expressed  as  the  sum  of  the  scattered  spherical  harmonics,  (Eq.(4)). 

It  is  interesting  to  reexpand  the  integrals  of  Eqs.(9)  -  (12)  in  terms  of  incoming 
spherical  harmonics  jn(k  R)P"1  (cos  (p).  In  order  to  do  this,  the  relation 

OO 

Jm(hr)e'^  =  i{n-m)P™(c°s4>)P'n('r/k)3n(kR),  (17) 

n~~m 

is  used.  Then  these  integrals  can  be  written  as 

°°  r°°  np2^D  R(n) 

r.  =  E  l  t7  1  +  w,),  s»P?h/k)Kn(-r/k)d, 

(18) 

00  r°°  np2i^h  —  R(n) 

r’  =  £i'+""2”'*<*s>  C<«»*>  I 

(19) 

°°  r°°  np2^^l~^)  _ l\ 

r3  =  E  r*"-*"3.V'W r(cos«  /  1  + 

v—rri  V”/ 

(20) 

Tt  =  E .:*+”-2";,(M)^‘(coS0)  r 

v—m 

(21) 

From  these  sums,  it  can  be  seen  that  an  outgoing  Spherical  Harmonic  function 
hn(kR)P™ (cos  (f>)  is  converted  into  a  set  of  incoming  Spherical  Harmonics  from  the 
reflections  off*  the  surrounding  boundaries.  After  truncating  the  set  of  outgoing  and 
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incoming  harmonics,  this  conversion  process,  at  order  m  can  be  expressed  by  a  con¬ 
version  matrix, 

vin{n)  =  Cmvout(n),  n  =  m, N  (22) 

where  vtn  is  the  vector  of  incoming  coefficients,  v^t  is  the  vector  of  scattering  coef¬ 
ficients  for  hn(kR)Py(cos(f))  and  vtnc  will  represent  the  original  incident  coefficients. 
The  (u,  n),  (u,  n  >  m)  element  of  the  conversion  matrix  Cm  is  computed  by  taking  the 
(v-m+1  )th  term  of  T\  +  T2  +  T3  +  T4  in  Eqs.(18)-(21)  which  represent  the  conversion 
of  the  outgoing  Spherical  Harmonic,  hn(kR)P™(cos  <f>)  . 

In  free  space,  the  scattering  process  for  the  harmonies  of  the  m th  azimuthal  order 
can  be  expressed  with  a  diagonal  scattering  matrix  D  whose  entries  are  the  scattering 
coefficients  aU)n  =  ra,  However,  because  of  the  additional  conversion  process 

described  above,  the  effective  scattering  from  the  sphere  is  given  by  a  sum  of  a  series 
of  scattering  and  reseat tering, 

v™*  =  {D  +  DCmD  +  DCmDCmD  +  (23) 

or 

v™1  =  ( D~l  -  Cmylvinc  (24) 


If  one  ignores  the  conversion  matrix  Crn  (the  single  scatter  approximation)  then 
Eq.(24)  reduces  to  vout  —  Dvtnc.  In  the  computation  of  the  scattering  function 
Fm(q)  of  Eq.(15)  one  can  ignore  the  reseattering  of  energy  onto  the  sphere  by  the 
surrounding  boundaries  and  simply  use  the  coefficients  an,n  =  m1...,oc  or  one  ean 
compute  the  effective  new  coefficients  using  Eq.(24). 


In  the  ease  that  the  incident  field  is  a  plane  wave  then  vinc(n)  can  be  specified  an¬ 
alytically.  However,  in  this  report,  the  incident  field  is  taken  to  come  from  a  point 
source  in  the  waveguide.  In  this  case,  the  m th  azimuthal  component  of  the  incident 
field  can  be  written  as 


00 

p'nc(r,  z ;  m)  -  ik  P”l(c os  0S)F"‘( cos  <p)hn(kRs)jn(kR)  +  i  / 

Jo 

tx{q)e^  +  T2(q)e +  f3{q)e~'yz  +  f4(q)e^z 


1  +  R{q)ei2^h 


qdq 


Jm{qrs)J,n{qr) 


(25) 


where  R,  and  (frs  are  the  source  coordinates  meaured  from  the  centre  of  the  sphere 
and  the  terms  T\,  T2,T3,  T4  are  given  by 


Tx(q) 

=  R{q)exp(i-f(D  4-  zrec)) 

(20) 

T2(h) 

=  -R{q)  exp(z7(D  +  2 h  -  zrec)) 

(27) 

T3(q) 

=  —  exp(z7(2/i  —  D  —  zrec)) 

(28) 

TM 

—  exp(f7(^rcc  d-  2 h  -^)) 

(29) 
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Using  the  relation  from  Eq.(17),  we  can  write  that  the  incident  term  for  P™( cos  <p)jn(kR) 
is  given  by 


Pnm  2 


f 


(fl(q)  +  f4(q))P?('y/k)  +  (f2(q)+f3(q))P^(-1/k) 


7(1  +  R{q)e2i^h) 


+  ikP™  (cos  <j>s)hn(kRs) 


J jn  ( qr, )  ([([([ 
(30) 


This  expression  when  combined  with  those  of  Eqs.(15)  and  Eqs.(9)  -  (12)  could  be 
expressed  as  double  integrals  and  expressions  similar  to  those  of  this  report  can  be 
found  expressed  in  Ref.9.  However,  in  this  report  the  scattered  field  is  computed  in 
the  manner  outlined  above.  First,  the  incident  coefficients  are  computed  using  the 
expression (s)  given  by  Eq.(30).  These  coefficients,  can  then  be  modified  by  Eq.(24), 
if  the  rescattering  effects  are  desired.  The  resulting  scattering  coefficients  are  used  in 
the  computation  of  Fm(h)  in  Eq.(15)  and  the  evaluation  of  the  scattering  integrals  of 
the  form  of  Eq.(14).  Once  the  scattered  fields  P^{r,z)  have  been  computed  for  all 
azimuthal  orders  m,  the  final  result  is  computed  from 

OO 

P(r ,  z,e)  =  ^2  (rnP,n(r,  z)  cos  mO  (31) 

m- 0 

where  in  our  examples  0  =  0. 


2.3  Scattering  from  a  sphere  in  a  waveguide  -  a  multipath 
expression 


The  multipath  expansion  expression  for  propagation  in  a  Pekeris  waveguide  is  well 
known  [12]  and  we  can  write 


P(r,  z) 
where 


J2L  plkDni  plkDn2 

+  EV>W-H  +  VjW-tj—  +  l-'j(n) 


pikDn3 


gik  Dn  4 


Aii  =  \/r2  +  (z  +  z~,  +  2 (n  -  l)/i)2  (33) 

D„  2  =  yfr2  +  (z  —  zs  —  2  n/i)2  (34) 

Dn 3  =  y/r2  +  (-z  -  zs  +  2 nh)2  (35) 

Ai4  =  \A2  +  (^  -  ~a  +  2n/i)2  (36) 


and  Do  is  the  distance  of  the  direct  arrival.  The  coefficients  are  given  bv 


Vi(n)  =  ir,(tfi„)(-l)"-1  (37) 

V2  (n)  =  (~R(0  2n))n  (38) 

V3  (n)  =  (  1)”-R(^3r»)”— 1  (39) 

V4(n)  =  (-P(M)"  (40) 
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Here  the  angles  f)\„, 02„,  etc.  denote  the  angle  that  the  ray  makes  with  bottom  seabed 
and  they  are  given  by 

9jn  =  Acos{-^~).  (41) 

LJjn 

and  R  denotes  the  seabed  reflection  coefficient  for  this  angle.  Wc  will  consider  N  rays 
of  each  of  the  4  ray  families  in  addition  to  the  direct  ray,  resulting  in  4iV  +  1  rays 
and  these  will  simply  be  denoted  with  an  index  k  =  1,  ...,47V  +  1  with  k  =  1  being 
the  direct  term. 


Using  this  approach,  the  rays  incident  upon  the  sphere  from  the  specified  source 
position  are  computed.  Then  the  same  geometric  set  of  rays  are  considered  the 
propagate  the  scattered  field  from  the  sphere  back  to  the  source.  This  approach  is 
correct  for  the  monostatic  case.  For  the  bistatic  case,  this  second  set  of  rays  would 
correspond  to  those  joining  the  sphere  centre  to  the  specified  receiver  position(s). 
Thus  wc  can  write  for  the  scattered  field, 


4jV  +  l  4jV  +  1 


P“(r,;)=  £  £  fWVW 


n\—  0  ri2— 0 


cxp{ik(Dnl  +  Pn2)) 
Dn  i  Dn  2 


S(0„2  -  0ni)-  (42) 


In  Eq.(42)  S((f>)  denotes  the  farfield  free  space  scattering  pattern  for  the  sphere  and 
due  to  the  symmetries  of  the  sphere,  this  function  depends  upon  only  the  one  angle, 
the  difference  between  the  incoming  and  outgoing  angles.  It  is  important  that  these 
angles  are  defined  in  such  a  way  that  the  difference  is  correct.  For  example,  for  the 
ray  scattered  upwards  from  the  sphere  and  then  reflecting  off  the  top  pressure  release 
surface  for  the  n3  family  of  Eq.(35),  the  angle  is  given  by 


=  Acos((2nh  —  z  —  zs)/Du 3)  (43) 

whereas  for  a  seabed  reflected  ray  of  Eq.(33),  the  angle  is  given  by 

=  7r  —  Tros((z  +  zs  +  2 (n  -  1  )h)/Dn\).  (44) 

The  far  field  scattering  pattern  S(cp)  is  given  by 


S(4>)  =  -if^aHX&^(-irPn(coB4>))  (45) 

n=0 

Although  the  multipath  approach  is  pleasing  intuitively,  it  is  not  clear  that  it  is 
mathematically  correct.  However,  as  is  done  in  [12],  in  the  wavenumber  integrals  in 
the  previous  expansion,  the  term  1/(1  -f  R(q)et2lh )  can  be  expanded  as 

1/(1  +  R(q)ei2^1)  =  1  -  R(q)et2lh  +  R2(q)e,^h  +  ...  (46) 
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Each  of  these  new  integrals  can  be  treated,  using  the  method  of  Stationary  Phase 
and  the  asymptotic  behaviour  of  Jm,  to  obtain,  in  the  end,  the  ray  expansion  just 
described  (note,  for  example,  in  Eq.(15)  that  we  wrote  the  factor  Frn(q)  in  terms  of  the 
scattering  pattern  for  the  sphere).  Alternatively,  one  could  use  the  same  approach  as 
Ingenito  [13].  Although  he  used  a  modal  representation,  he  interpreted  the  scattering 
in  terms  of  the  modes  having  equivalent  ray  directions  and  being  a  combination  of 
plane  waves. 

In  the  monostatic  case,  Eq.(42)  can  be  rewritten  as 


4N+1  4N+1 


exp  (ikjDrg  +  Pn2)) 
Dn\Dn  2 


S{<j)n2  ~  <i>nl),  (47) 


«1=0  ri2=ni 


where  tu(ni,n2)  =  1  for  rt\  =  n2  and  w(ni,n2)  =  2  otherwise.  Although,  we  will 
not  consider  a  moving  source  in  this  report,  the  form  of  Eq.(47)  is  very  amenable  to 
considering  a  doppler  shift  for  each  combination  of  ray  path. 
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3  NUMERICAL  EXAMPLES 


We  start  by  considering  scattering  from  an  elastic-shelled  and  a  rigid  sphere  in  free 
space.  The  sphere  is  0.25  m  in  radius.  The  shell  is  G  mm  thick  and  is  steel.  The 
interior  of  the  sphere  is  taken  to  be  a  vacuum.  The  scattering  coefficients  an  are 
eomputed  for  8001  frequencies  at  an  increment  of  10  Hz  for  both  the  shelled  and  rigid 
case.  The  rigid  ease  eorresponds  to  the  normal  derivative  of  the  pressure  vanishing  at 
the  surface  of  the  sphere.  In  this  case,  the  scat  tering  coefficients  are  given  analytically 
by 


wi 


(48) 


The  coefficients  in  the  ease  of  an  elastie-shelled  sphere  are  eomputed  using  an  existing 
FORTRAN  eode.  For  the  case  of  the  sphere  in  a  waveguide,  these  coefficients  are 
read  in  as  a  file  for  the  MATLAB  code  which  then  uses  them  in  the  waveguide 
computations  using  either  the  wavenumber  integral  approach  or  the  multipath-ray 
approach.  In  Fig.  1  the  spectrum  of  the  backscattered  field  at  a  range  of  100  m 
is  shown  for  the  shelled  sphere  in  blue  and  for  the  rigid  sphere  in  red.  There  are 
strong  resonance  features  for  the  shelled  sphere.  In  Fig.  2  we  show  the  computed 
backscattered  timeseries  for  an  incident  pulse  of  the  form 


S(f)  =  exp (-(/  -  fc)2/(2(B\V/2)2)),  f  >  0  (49) 


where  fc  —  (/ 1  +  fen<i)/ 2  and  BW  =  (fend  ~  /i)-  For  all  the  pulses  f\  —  100// z 
with  /r„d  =  4100  Hz,  8100Hz,  and  16100  Hz.  The  left  hand  column  shows  the  time 
series  for  the  elastic  shelled  sphere  and  the  right  hand  column  shows  the  response  for 
the  rigid  sphere.  It  is  elear  that  there  is  significant  return  for  the  shelled  sphere  for 
several  milliseconds  after  the  specular  return,  particularly  for  the  lowest  frequency 
pulse  centred  at  about  2  kHz.  In  Fig.  3,  a  zoom  of  the  region  near  the  specular 
reflection  is  shown  for  the  elastic  and  rigid  spheres  for  the  [100  16100]  Hz  pulse.  The 
specular  reflection  for  the  shelled  sphere  has  the  opposite  polarity  than  that  of  the 
rigid  sphere.  This  is  because  the  shell  is  still  relatively  thin  at  this  frequency  and  the 
interior  of  the  sphere  is  a  vacuum.  The  rigid  sphere  has  a  small  secondary  arrival 
after  the  speeular  reflection  -  this  is  the  Franz  wave  arrival  [14]. 

The  scattering  coefficients  an(uj)  for  the  elastic  shelled  sphere  are  now  used  in  con¬ 
junction  with  Eqs.(15)  and  (16)  to  compute  the  scattered  field  using  the  wavenumber 
integral  approach  and  also  with  the  ray-based  approach.  The  seabed  has  a  sound 
speed  of  1700  m/s  and  a  density  of  p  =  1.5  g/cm3.  The  depth  of  the  oeean  waveg¬ 
uide  is  30  m.  The  source  is  located  at  a  depth  of  15  m.  The  sphere  is  located  at  a 
horizontal  range  of  100  m  and  its  depth  is  varied  between  0.3  to  29.7  m.  The  scat¬ 
tered  field  is  computed  for  2001  frequencies  in  the  range  [100  4100]  Hz.  In  Figs.  4  and 
5,  we  show  a  comparison  of  the  computed  backscattered  spectra  (multiplied  by  104 
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Figure  1:  The  backscattered  spectrum  (free  space)  for  an  elastic-shelled  sphere  (red), 
0.25  m  in  radius  and  a  rigid  sphere  (blue) 


to  compensate  for  range)  for  sphere  centres  depths  ranging  from  0.25  in  to  29.75  m. 
The  depth  of  0.25  m  means  that  the  sphere  is  just  touching  the  upper  pressure  re¬ 
lease  surface  and  the  depth  of  29.75  m  signifies  the  sphere  touching  the  seabed.  In 
Figs.  6  and  7,  the  resulting  time  series,  using  the  source  function  of  Eq.(49)  with 
a  bandwidth  of  4  kHz,  are  shown.  As  can  be  seen,  the  inclusion  of  the  rescattering 
effects  has  a  significant  effect  for  the  sphere  close  to  the  upper  surface.  This  effect 
dies  away  rapidly  as  the  sphere  centre  goes  more  than  about  40  cm  from  the  surface. 
For  this  particular  example,  the  rescattering  terms  do  not  have  a  significant  effect  for 
the  sphere  close  to  the  seabed.  However,  it  is  not  clear  that  this  will  always  be  the 
case.  In  fact,  since  the  bottom  seabed  can  be  made  to  approximate  either  a  pressure 
release  or  rigid  surface,  this  will  definitely  not  always  be  the  case. 

In  the  following  computations,  we  will  ignore  the  rescattering  terms.  However,  it 
should  be  noted  that  for  the  sphere  within  small  distances  of  the  free  surface  (and 
perhaps,  the  seabed)  these  results  may  not  be  accurate.  It  is  possible  to  include 
these  effects  in  the  ray-based  model  by  including  these  terms  and  computing  a  two- 
dimensional  scattering  function  5(</w,  </>tn)  (the  scattering  function  would  no  longer 
simply  depend  upon  the  difference  of  angles)  .  In  the  case  of  varying  target  position 
this  would  be  time-consuming,  as  the  rescattering  coefficients  would  need  to  be  re¬ 
computed  for  each  new  target  position.  In  the  case  when  the  target  position  is  fixed 
and  the  receiver  or  source  position  varied,  this  might  be  a  feasible  approach. 

In  Fig.  8  we  show  a  comparison  between  the  results  obtained  from  the  wavenumber 
method  (no  rescattering)  and  the  multipath  expansion  method.  Here,  at  100  Hz, 
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Figure  2:  The  computed  backscattered  time  series  for  a  Gaussian  weighted  spectrum 
-  elastic  response  shown  in  (a)-(c)  for  frequencies  in  [100  4100]Hz,  [100  8100]  Hz  and 
[100  16100]  Hz.  The  corresponding  rigid  results  are  shown  in  (d)-(f) 


Time(ms) 


Figure  3:  A  detailed  look  at  the  region  of  the  specular  reflection  for  Figures  (c)  and 
(f)  of  Fig.3  for  the  elastic  shelled  and  rigid  spheres.  The  frequency  range  is  [100 
16100}  Hz. 
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Figure  4:  A  comparison  of  the  spectra  as  a  function  of  sphere  depth  with  rescattering 
terms  (blue)  and  without  rescattering  (red) 


Figured:  A  comparison  of  the  spectra  as  a  function  of  sphere  depth  with  rescattering 
terms  (blue)  and  without  rescattering  (red) 
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Figure  6:  A  comparison  of  the  time  series  as  a  function  of  sphere  depth  with  rescat¬ 
tering  terms  (blue)  and  without  rescattering  (red) 
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Figure  7:  A  comparison  of  the  time  series  as  a  function  of  sphere  depth  with  rescat¬ 
tering  terms  (blue)  and  without  reseat teriug  (red) 
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Figure  8:  The  scattered  held  amplitude  (multiplied  by  IQ4)  as  computed  by  wavenum¬ 
ber  integral  approach  (blue)  and  the  multipath  approach  (red)  as  a  function  of  sphere 
depth.  The  source/ receiver  depth  is  15  m  and  the  range  to  the  sphere  is  100  m. 


there  is  a  significant  disagreement  between  the  results.  This  is  not  surprising.  The 
multipath  expansion  approach  is  a  high  frequency  approach.  For  example,  the  be¬ 
haviour  of  the  wavenumber  integral  near  the  critical  point  is  not  accounted  for  and 
effects  such  as  beam  displacement[15]  are  ignored.  The  multipath  curve  shows  2  slight 
“kinks”  near  8.3  and  21.7  m.  These  correspond  to  depths  for  which  there  are  2  bot¬ 
tom  reflections  with  an  angle  very  close  to  the  critical  angle  for  this  seabed.  However, 
by  1000  Hz,  the  2  results  are  nearly  identical.  In  Fig.  9  we  show  a  comparison  of  the 
2  curves  for  100  -  400  Hz.  In  this  case,  there  is  still  disagreement  between  the  curves 
at  200-400  Hz,  but  they  are  much  closer  than  at  100  Hz. 

In  Fig.  10  we  show  a  two-dimensional  plot  of  the  received  spectra,  computed  using 
the  multipath  expansion  method  as  a  function  of  frequency  and  sphere  depth  and 
in  Fig.  11  we  show  the  resulting  time  series  computed  using  Fourier  synthesis.  In 
both  plots  one  can  see  the  significant  effect  of  the  sphere  depth  on  the  character  of 
the  received  signal  and  also  the  semi-periodic  appearance.  In  Fig.  12  we  show  in  (a) 
the  direct  echo  pulse,  in  (b)  the  time  series  for  a  sphere  at  0.3  m  depth  (c)  for  the 
sphere  at  15  m  depth  and  (d)  the  sphere  at  29.7  m  depth  (i.e.,  just  above  the  seabed). 
Once  again,  the  large  variation  of  the  echo  with  sphere  depth  is  seen  as  well  as  the 
significant  change  from  the  echo  which  would  be  received  in  free  space. 

We  now  show  the  spectral  plot  and  time  series  for  the  case  of  the  sphere  at  600  m 
range.  It  is  interesting  to  consider  how  many  multipath  terms  are  required  to  obtain 
an  accurate  pressure  computation.  Let  us  suppose  that  6C  denotes  the  critical  angle 
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Figure  9:  The  scattered  Held  amplitude  (multiplied  by  104j  as  computed  by  wavenum¬ 
ber  integral  approach  (blue)  and  the  multipath  approach  (red)  as  a  function  of  sphere 
depth .  The  source/ receiver  depth  is  15  m  and  the  range  to  the  sphere  is  100  m. 
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Figure  10:  The  scattered  field  amplitude  (scaled  by  104)  as  function  of  frequency  and 
sphere  depth 
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Figure  11:  The  scattered  held  time  series  (amplitude)  as  a  function  of  sphere  depth. 
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Figure  12:  Representative  time  series  for  (a)  the  free  space  echo  (b)  the  sphere  at  10 
depth  (c)  the  sphere  at  15  m  depth  and  (d)  the  sphere  at  29.7  m  depth 
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Figure  13:  The  scattered  field  amplitude  (scaled  by  C  *  104J  as  function  of  frequency 
and  sphere  depth  for  a  range  of  600  m . 


for  the  seabed,  measured  off  the  horizontal.  Then  the  cycle  distance  corresponding 
to  that  angle  is 

rc  =  2  h/tan(9r)  (50) 

For  our  example,  the  cycle  distance  is  about  112  m.  For  the  multipaths  ,  which 
correspond  to  steeper  angle  than  this,  they  will  quickly  die  out  due  to  the  power 
of  the  reflection  coefficient,  which  is  less  than  one  in  this  case.  Thus,  for  a  range 
of  600  m,  approximately  6  terms  should  suffice.  To  be  on  the  safe  side,  we  will 
consider  8.  In  the  previous  example,  at  100  m  range,  we  simply  used  the  frequency 
spacing  of  10  Hz  at  which  the  free  space  scattering  coefficients  were  computed.  At 
longer  ranges,  this  frequency  spacing  is  not  sufficiently  fine,  due  to  the  larger  travel 
time  differences  of  the  paths.  In  Figs.  13-14,  we  show  the  resulting  spectral  and 
time  pulse  computations  for  a  range  of  600  m.  The  time  series  have  been  scaled 
by  a  factor  of  6  which  makes  their  amplitudes  approximately  equal  to  those  for  the 
100  m  case.  In  this  approach  we  simply  took  a  brute  force  approach  of  computing 
the  scattering  coefficients  ak  at  a  spacing  of  1  Hz.  However,  in  order  to  speed  the 
computations  for  subsequent  computations,  we  have  rearranged  the  original  code 
(which  was  vectorized  with  respect  to  the  sphere  position)  so  that  for  a  fixed  sphere 
depth  and  a  particular  multipath  arrival/scatter  type,  all  the  frequency  computations 
are  done  in  a  vectorized  form.  Many  of  the  terms  are  frequency  independent  and  need 
be  computed  only  once. 

Finally,  we  consider  a  source  with  a  bandwidth  of  80  KHz.  The  spectrum  of  the  direct 
return  is  that  of  Fig.  1  (weighted  by  the  spectrum  of  the  source).  The  range  is  100  m 
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Figure  14:  Representative  time  scries  (scaled  by  a  factor  of  6)  for  (a)  the  free  space 
echo  (b)  the  sphere  at  10  m  depth  (c)  the  sphere  at  15  m  depth  and  (d)  the  sphere 
at  29.7  m  depth  for  a  range  of  600  m. 

but  now  the  source/receiver  is  2  m  above  the  seabed.  In  Fig.  15  we  show  the  resulting 
spectra  as  a  function  of  frequency  and  sphere  depth  for  the  first  [100  4100]  Hz, 
corresponding  to  the  first  example  we  considered,  but  with  a  different  source/receiver 
depth.  This  is  a  “zoom  ”  of  this  region  contained  within  the  much  larger  [100  80100] 
Hz  interval  shown  in  Fig.  17.  As  can  be  seen,  the  interference  patterns  seen  in  this  case 
are  quite  different  then  the  source/receiver  being  a  midwater  depth.  Some  selected 
time  series  are  shown  in  Fig.  1G.  For  the  case  of  the  sphere  on  the  seabed,  it  is 
predicted  that  there  will  be  virtually  no  direct  arrival  energy  due  to  the  constructive 
interference  of  the  seabed  reflected  energy.  The  resulting  backscattered  spectra  as  a 
function  of  sphere  depth  and  frequency  are  shown  in  Fig.  17  for  the  entire  80  kHz 
bandwith.  The  resulting  time  series,  for  a  Gaussian  pulse  centred  at  40.1  KHz  and 
a  bandwidth  of  20  kHz,  are  shown  in  Fig.  18.  For  these  higher  frequencies  and 
larger  bandwith,  the  individual  arrivals  seem  better  resolved.  However,  there  is  still 
confusion  because  each  multipath  arrival  consists  of  its  own  sequence  of  arrivals,  due 
to  the  elastic/structural  scattering  of  the  sphere.  In  Fig.  19  a  zoom  of  the  time  series 
near  the  specular  reflection  are  shown  for  (a)  just  the  direct  paths  and  (b)  for  all 
paths,  for  the  case  of  the  sphere  centre  0.3  mn  above  the  seabed.  It  can  be  seen  in 
(b)  that  the  first  echo  appears  to  be  different  due  to  the  coherent  combination  of 
the  direct  and  the  seabed  reflected  paths.  This  alteration  of  the  specular  echo  was 
observed  with  real  data  and  discussed  in  Ref.  5. 
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Figure  15:  The  scattered  field  amplitude  (scaled  by  101)  as  function  of  frequency  and 
sphere  depth  for  a  range  of  100  m  and  source/receiver  at  2  m  above  the  seabed. 


Time(sec) 

Figure  16:  Representative  time  series  for  (a)  the  free  space  echo  (b)  the  sphere  at 
10  in  depth  (c)  the  sphere  at  15  m  depth  and  (d)  the  sphere  at  29.7  in  depth  for  a 
range  of  100  m.  The  source  function  was  a  Gaussian  pulse  centred  in  [100  4100]  Hz. 
The  source/ receiver  are  at  2  m  above  the  seabed. 
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Figure  17:  The  scattered  field  amplitude  (scaled  by  104)  as  function  of  frequency  and 
sphere  depth  for  a  range  of  100  m  and  source/ receiver  at  2  in  above  the  seabed. 


Time(sec) 

Figure  18:  Representative  time  series  for  (a)  the  free  space  echo  (h)  the  sphere  at 
10  m  depth  (c)  the  sphere  at  15  in  depth  and  (d)  the  sphere  at  29.7  m  depth  for  a 
range  of  100  m.  The  source  function  was  a  Gaussian  pulse  centred  at  40100  Hz  with 
a  standard  deviation  of  20  kHz.  The  source/receiver  are  at  2  m  above  the  seabed 
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Figure  19:  .  A  zoom  of  the  echo  for  (a)  direct  paths  only  (b)  all  paths  for  the  region 
of  the  specular  reflection.  The  source  function  was  a  Gaussian  pulse  centred  at  40100 
Hz  with  a  standard  deviation  of  20  kHz.  The  source/receiver  are  at  2  m  above  the 
seabed 


4  SUMMARY 


We  have  shown  in  this  paper  how  the  scattering  by  a  spherical  elastic  shell  can 
be  accurately  modelled  within  a  simple  Pekeris  waveguide.  It  was  found  that  the 
rescattering  terms  of  the  sphere  with  the  surrounding  boundaries  could  be  significant 
for  the  sphere  very  close  to  a  boundary  (in  our  example,  the  upper  surface).  It 
was  found  that  for  frequencies  greater  than  about  500  Hz,  that  a  multipath/iinage 
solution  could  be  used  to  model  the  propagation  and  scattering.  The  propagation 
effects  within  the  waveguide  and  the  sphere’s  depth  in  the  waveguide  alter  the  received 
echo  significantly.  In  future  work,  we  will  use  this  model  to  examine  the  resulting 
possibilities  for  target  classification. 
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combinations  of  incident  and  backscattered  multipath  arrivals.  For  a  very  shallow  target  or  one  close 
to  the  seabed,  it  may  be  difficult  to  isolate  some  of  these  multipath  arrivals.  It  is  expected  that  these 
waveguide  effects  could  have  a  significant  effect  upon  some  classifiers.  In  this  paper,  we  discuss  the 
modelling  of  scattering  of  sound  from  an  elastic  sphere  (shelled,  solid,  water-filled,  etc)  which  is  in  a 
waveguide  with  an  upper  surface  and  a  lower  seabed.  An  exact  propagation  solution  is  implemented 
and  used  to  benchmark  a  ray-expansion  solution.  The  scattered  signals  for  a  variety  of  different  sonar 
pulses  and  different  positions  of  a  spherical  target  in  the  waveguide  are  shown. 
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